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Abstract 

The Pariser-Parr-Pople Hamiltonian is used to calculate and identify the nature of 
the low-lying vertical transition energies of polydiacetylene. The model is solved using 
the density matrix renormalisation group method for a fixed acetylenic geometry for 
chains of up to 102 atoms. The non-linear optical properties of polydiacetylene are 
considered, which are determined by the third-order susceptibility. The experimental 
l_B,j data of Giesa and Schultz are used as the geometric model for the calculation. 
For short chains, the calculated E{lBu) agrees with the experimental value, within 
solvation effects (~ 0.3 eV). The charge gap is used to characterise bound and unbound 
states. The n5„ is above the charge gap and hence a continuum state; the IBu, 2Ag 
and mAg are not and hence are bound excitons. For large chain lengths, the nBu 
tends towards the charge gap as expected, strongly suggesting that the nBu is the 
conduction band edge. The conduction band edge for PDA is agreed in the literature 
to be ~ 3.0 eV. Accounting for the strong polarisation effects of the medium and 
polaron formation gives our calculated EooinBu) ~ 3.6 eV, with an exciton binding 
energy of ~ 1.0 eV. The 2Ag state is found to be above the IBu, which does not 
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agree with relaxed transition experimental data. However, this could be resolved by 
including explicit lattice relaxation in the Pariser-Parr-Pople-Peierls model. Particle- 
hole separation data further suggest that the IBu, 2Ag and mAg are bound excitons, 
and that the nBu is an unbound exciton. 



I INTRODUCTION 

Polymers, and other molecular materials, that exhibit non-linear optical properties and elec- 
troluminescence have attracted much interest amongst theorists and experimentalists, owing 
to their possible use in organic technology [Q-^. Poly(diacetylene)s (PDAs) with the gen- 
eral formula shown in Fig. 1. (a) are the vr-conjugated organic polymers studied here: they 
exhibit non-linear optical properties, near-perfect crystal structure and doping-dependent 
transport properties. These properties make them ideally suited for use in the theorist's 
calculation. A full characterisation of PDA's electronic properties would mean a better un- 
derstanding of TT-conjugated polymers in general. In PDA with a very low polymer content 
{xp < 10~^ in weight), interchain interaction is minimal and thus it can be considered an 
ideal one- dimensional model system. 

There have been many calculations performed on PDA with varying degrees of suc- 
cess. Parry conducted a SCF calculation and found an energy gap that was too low 
Boudreaux's SCF-Xa calculation gives an energy gap of 2.21 eV, which agrees well with the 
long-chain experimental extrapolations for the energy gap of 2.11 eV and 2.25 eV (yellow 



phase) [3-|[T0|. A further SCF-Xa calculation by Boudreaux [|Tl|| accounts for charge density 



waves by the same geometric Hiickel model of Whangbo [T^. The ab initio crystal orbital 



calculations of Kertesz give energy gaps that are too high [jT3[-|jTj]. Predicted bond lengths 



show impressive agreement with x-ray data in Karpfen's ab initio SCF calculation |15]. He 



performed a total energy minimisation to calculate bond lengths and shows that the energy 
difference between butatriene and acetylene structure per unit cell is 0.52 eV. This agrees 



with previous calculations and experiment |T^, |T6[. A similar calculation yielding similar 
results was performed by Bredas |[T^ . 

Wilson, Cojan, Cade and Movaghar ||l8[-p0[ performed Hiickel band structure calcula- 
tions for the infinitely long PDA chain. Cajon and Wilson did not include bond alteration 
in their calculation, and thus modelled a compound with one distinct chemical bond. The 
energy gap was found to be too low. Cade and Movaghar used two degrees of freedom in 
their calculation bond alteration parameter. 



Extended Hiickel calculations have been performed by Parry and Whangbo [O, pT|. The 



work of Parry is restricted to one distinct chemical structure and, as with other parameter- 
isations, yields a value for the energy gap that is too low. The model employed by Wangbo 



has a hopping term of the form: 



t = toil±S) 



(1) 



with < 6 < 1. This produces results that are, again, quantitatively too low; however, it does 
generate an energy difference of 0.48 eV between the butatriene and acetylene structures, in 
excellent agreement with other calculations and experiment |T^-|T^. In summary, therefore, 
all these models are seen to give good qualitative predictions for bond alteration, but the 



data are seen to be consistently red-shifted owing to their single particle picture. 

The importance of screening in PDA is still considered controversial to many authors, 
because it determines the role of electron-electron correlations. A simple way to account for 
screening is to adjust the single chain electron-electron interaction potential V{r), a method 
suggested by many authors. For example, Abe lowered the Hubbard parameter U from 
11 to 5 eV and used a large long-range dielectric constant in a Ohno potential p2|-p3 



However, this again causes red-shifted energy gaps and results in unphysical atomic cohesion 
energies. We emphasise that there is clear and indisputable evidence that electron-electron 
interactions are of paramount importance in the electronic structure calculations of PDA. 
First, in polymeric crystal form, the linear absorption spectrum is see to be symmetric and 
peaks at ~ 1.8-1.9 eV. However, the onset of photoconduction occurs not at this value 
but at ~ 2.3-3.0 eV (corresponding to uncorrelated vr-vr* transitions) p4|. In addition, 
comparison of the linear absorption with the EA spectrum shows an EA feature at ~ 2.0 
eV corresponding to the first derivative of the absorption: this is attributed to the Stark 
shift of the exciton. There is also an oscillating Franz-Keldysh band-edge structure at ~ 
2.5-3.0 eV, with a lineshape that deviates from the first derivative of the absorption. The 
oscillations are due to interband transitions and coincide with the onset of photoconduction. 



These data independently suggest that the exciton binding energy is ~ 0.5 eV [^-[^. We 
note that these values show a great deal of variability owing to the variety of phases and 
the disorder present in PDA. Hence, despite the apparent success of the Su-Schrieffer-Heeger 



model [^-[^ neutral excitons with large binding energies are the low-lying excitations in 
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PDA p4|, This is borne out by the Pariser-Parr-Pople calculation presented in 

this work, which incorporates electron-electron interactions. 



II MODELLING POLYDI ACETYLENE 
A Electronic Structure 

The atomic orbital model used in the calculation is described as follows. Carbon has the 
atomic configuration ls^2s^2p^. Since carbon has four valence electrons, there is an alter- 
nating structure of two sp^ and two sp hybridised carbon atoms that can account for the 
idealised PDA chemical composition. (See Fig. 1.) 

Since PDA is spatially centrosymmetric, and thus shows C2 symmetry, the wavefunctions 
possess mirror plane and centro-inversion symmetries. The group notation for mirror plane 
symmetries are A for the symmetric and B for the antisymmetric case. Inversion symmetries 
are labelled g for symmetric and u for antisymmetric. The ground state is therefore labelled 
lAg and the first optically active dipole from the ground state has to be the state. The 
wavefunctions are either even (^4^) or odd {By) under an inversion. 

B Geometric Structure 

Although there is a myriad of literature on the optical properties of PDA, there is still 
little on their structural and conformational properties, which prevents extensive theoretical 
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investigation. The work of Giesa and Schultz following that of Wudl and Bitler |32 
describes the synthesis and thorough characterisation of the series of alternating all-trans- 
polyenynes without substituents at the vinylic bonds that are studied in this paper. 

To maintain consistency with the literature, the unsubstituted model compounds are 
named as follows: 



Cn[N] (2) 

where the number of unsaturated carbon atoms n specifies each compound, and the number 
of formal monomer units N are added in the square brackets for clarity. The relation between 
n and N is. 



n = 2(2iV + l). (3) 

These structures are shown in Fig. 2. In contrast with the simplest conjugated polymer 
represented by trans-p olyacetylene, the carbon backbone of PDA contains two, additional, 
localized vr electrons in every unit cell. Thus, there are two possibilities of bond alternation 
in this system that lead to non-equivalent structures with non-degenerate ground states, as 
shown in Fig. 1. (b) and (c). Work on the molecular geometries of Ci4[3] and C22[5] by X-ray 
structure analysis has determined bond lengths and angles that represent a typical acetylene 
structure (see TABLE I.) and no evidence for a hutatriene form is found. In addition, Giesa 
shows that substitution does not alter the polymer structure significantly. Increasing the 
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size of the polymer does elongate double and triple bonds, while reducing the single bond; 
however, this is not pronounced. Hence, overall, these findings justify the use of the model 
acetylene structure employed in this paper, because an increase in conjugation length does 
not show a substantial effect on bond lengths and angles. Hence, owing to PDA's simple 
electronic and geometric ground-state structure, the Pariser-Parr-Pople Hamiltonian can 
easily describe it. 



Ill THE PARISER-PARR-POPLE HAMILTONIAN 

The Pariser-Parr-Pople Hamiltonian is: 

Hppp = - tiji4aCja + c]^Ci„) + U Y^i^H -hi^i-h 

<ij>a i 

+ E^*j(^- -1) (4) 

where the operator c\^{cia) creates (annihilates) a electron of spin a at site i (the zth 
carbon atom), tij (>0) is the transfer integral between the nearest neighbour atomic orbitalsl. 
The first term in the Hamiltonian allows the electrons to "hop" from one atom to another, 
and represents the kinetic energy gained from delocalising an electron from its atomic site. 
nia- = CiaC-^ is the number density operator, and < ij > denotes nearest neighbours. U 
is the Pariser-Parr-Pople Coulomb repulsion between two electrons occupying the same 2pz 



orbital and U = 10.06 eV [^. Vij is the long-range Coulomb repulsion and, in this work. 



is the Ohno function, Vij = U/ + Prfj where (3 = ([7/14.397)^ and rjj is the interatomic 
distances measured in A ||37l|. 



The general expression for the resonance integrals is given by: 

tij = to + a{rij - ro) (5) 

where a is the coupling constant of electron-phonon interactions, rij is the length of the bond 
between carbon atoms i and j, and tq is the average bond length of the reference system. 

The parameters in expression are derived using polyacetylene as a reference system. 
First, it is pertinent to consider the gradient of the resonance integral with respect to bond 
length. It has been shown that the difference 5t of the resonance integral is 0.2162 eV and the 



corresponding average change in bond lengths 61 is 0.0471 |3^. Thus, the electron-phonon 
coupling constant is: 

5t 0.2162 
"=57 = 0:0471 ^^ 

A full list of Hamiltonian parameters derived from polyacetylene used in this work is shown 
in TABLE II., and a schematic of the unit cell is found in Fig. 3. 

The different resonance integrals are determined by the specific structures of the chain 
with lattice constant 4a. Perturbations of the equidistant carbon chain lead to the energy 
gaps Eg and The ground state corresponds to a fully filled valence band and unoccupied 
conduction band. The IBu is reached from the ground state by exciting one electron from 
the highest occupied molecular orbital (HOMO) to the lowest unoccupied molecular orbital 
(LUMO) (see Fig. 4.). 

Although this band picture is useful in identifying important states, it has been shown 
that interactions do play an important part in the physics of PDA, as mentioned earlier. 

8 



By including the on-site and long range terms explicitly in the Hamiltonian electron- 
electron interactions are taken into account, giving an accurate treatment of the electronic 
structure of PDA. Turning on the interactions, U and Vij, changes the relative position of 
electronic states and can even cause the crossing of states. This non-interacting picture 
is therefore modified and an excitonic picture formulated. An electron and hole can bind 
together by their Coulomb field to produce an exciton, or bound electron-hole pair. Since 
the conduction band signifies the non- interacting limit, excitons are energetically situated 
below the conduction band. 



IV RESULTS AND DISCUSSION 



A Essential States 



Non-linear optical experimental data can be used to characterise the low-lying states of vr- 
conjugated polymers. The relative ordering of these states can elucidate much of the physics 
of such systems I^O], ^^j-]^. Non-linear optical processes in polymers with C2 symmetry are 



determined by the third-order susceptibility, ^^i — ^2 — ^3; ^1, ^2, 1^3) , which is calculated 
by the sum over all available states. This makes calculating x^^^ in principle difficult, owing to 
the large number of paths included. However, it is found that there are certain essential states 
that are important when considering the x^^^ oi 7r-conjugated polymers, as only the essential 
states make a significant contribution to it. The four essential states are the ground state, 
the first odd-parity exciton state, IBu, the charge transfer state, mAg, and the nBu (see Fig. 



5.). Another important state is the 2Ag: this is a two-photon state, whose position relative 
to the IBu determines whether a polymer exhibits photoluminescence. The importance of 
these states is the result of the strong and dominant dipole couplings amongst them in 

^^^f = ifim (7) 

resulting in their substantial contribution to non-linear optical spectroscopy. Hence, an 
important first step in identifying the non-linear properties of PDA is calculating the dipole 
moments of transitions between states. As the DMRG method is a robust and accurate 
method of finding dipole moments the essential states can be found. These are given in 
TABLE III. However, they are only calculated for polymers of up to 26 sites. This is 
because the important optical states become interlaced with other states of a spin-density- 
wave character (i.e., those related to the 2Ag). Since we can only target approximately 10 
states in each symmetry sector, the mAg and the nBu states soon become impossible to 
track. 

An analysis of the essential states and the 2Ag forms the remainder of this discussion. 
They are paramount when considering the physics of vr-conjugated polymers, as mentioned 
above, and are compared with experimental data. In Section |IV| (0) the vertical excitation 
energies of the essential states and the 2Ag are analysed. Section |^ (0) contains particle- 



hole separation data, which is used to confirm predictions made in Section |^ dBj) and help 
identify the nature of the IBu, 2Ag, nB^ and mAg states. 
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B Vertical Excitation Energies 



The DMRG method, used by Barford and Bursill following the work of White [Q, was 
implemented on the model system to find the vertical low-lying energy eigenvalues of PDA. 
The low-lying excitations are bound particle-hole pairs, which behave as composite particles. 
These particles delocalise along the polymer backbone as effective single particles, and thus 
their dispersion should be characteristic of single particles. In PDA the excitons effectively 
tunnel between double and triple bond dimers. In an effective particle model, this would 
be modelled as a linear chain, with two 'sites' per unit cell. The dispersion should scale as 
as these are states within the exciton band. In contrast, the band gap is associated 
with free, single particle transitions across the HOMO-LUMO gap. An examination of 



the dispersion relation of Lennard- Jones [P] indicates that this energy should scale as 1/n 
in the asymptotic limit. However, neither a solely 1/n nor 1/n^ fit for the free particle- 
hole and exciton dispersions, respectively, is accurate for small chains, owing to higher-order 
corrections, and the fit required cannot be used to characterise states. In this work, therefore, 
a polynomial fit is preferred for extrapolating the data. 

The calculated vertical excitation energies for the IB^, 2Ag, mAg, nB^ and charge gap, 
A, are plotted in Fig. 6. against 1/n. In addition, the polynomial extrapolations for the 
long-chain limit are shown in TABLE IV. for some of the states. The charge gap is a useful 
criterion for characterising excitonic and unbound states, and is given as follows, 

A = E(n + 1) + ^(n - 1) - 2E(n) (8) 
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Here E{n) is the ground-state energy of the n electron system. The charge gap signifies the 
lowest energy excitation of an electron from the valence to conduction band. Thus states 
above A are unbound excitons; those below it are bound. Of the states plotted in Fig. 6. 
only the nB^ is positioned above the charge gap, suggesting that the ni?„ is a continuum 
state. Conversely, the IB^, 2Ag and mAg, are below it and are thus bound excitons. This is 
shown to be true in Section |IV| (|C|) by examining the particle-hole separation of these states. 



The experimental data obtained by Giesa |TD| are also plotted in Fig. 6. A polyno- 
mial extrapolation of the experimental data yields i?oo ~ 2.5 eV; however, this is considered 
unreliable as there are so few points used in the fit. For short chains the calculated and 
experimental IBu energies are remarkably within solvation (polarisation) effects of the ex- 
perimental medium (~ 0.3 eV). This solvation value is derived from the work of Yaron, who 



predicts 0.3 eV for the IB^ [23|. In addition, Barford and Bursill's calculation on polyenes 



found E(lBy) to be 0.3 eV above the experimental value. Their recent calculations on 



PPP and PPV yield similar results [gSf. Our calculated E^{lBu) for PDA is ~ 3.0 eV, and, 
hence, correcting for solvation effects gives our calculated Eoo{lBu) ~ 2.7 eV. 

The energy of the nB^ is seen to tend to that of the charge gap for n — cxd, as seen in Fig. 
6. Hence, the extrapolated long chain nBy_ energy is found to be ~ 5.7 eV, and this state 
is predicted to be the conduction band threshold. However, the onset of photoconduction 
(continuum limit) is ~ 3.0 eV. The large polarisation effects of the surrounding medium can 
correct for this: these are estimated by Yaron to be ~ 1.5 eV [^]. In addition, polaron 
formation means a further ~ 2 x 0.3 eV can be subtracted from 5.7 eV [^. Including all 
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these effects gives Eoo{nBu) ~ 3.6 eV, which is in reasonable agreement with the experimental 
conduction band threshold quoted earlier. Therefore, the exciton binding energy is Ei, = 
Eoo(nBy) — Eoo(lBu), and from our calculation is found to be 3.6 — 2.7 ~ 0.9 eV. This is in 



reasonable agreement with binding energies found in the literature |25]-[27]. 

Having identified the important non-linear optical states in PDA from the dipole mo- 
ments, and calculated the corresponding excitation energies, an energy level diagram can be 
drawn. This is shown for 26 atoms in Fig. 7. As expected there is a large lAg IB^ dipole 
moment, corresponding to the one-photon absorption, while the band threshold {nBy) has a 
weak coupling to the ground state. It is seen that the first excited Ag state is above the IB^, 
which would suggest that PDA shows photoluminescence. (There is a direct lAg IB^ 
transition.) However, experiments have shown that PDA is not electro- luminescent: mea- 
surements of the frequency dependence of third harmonic generation in Langmuir-Blodgett 
films of PDA indicate that an Ag singlet state is ~ 0.11 eV below the lowest energy B^ state. 
In addition, the recent work of Kohler and Schilke has shown an Ag symmetry state ~ 
0.21 eV lower than the lowest B^ symmetry state. Our model calculates vertical transitions 
only and is seen to agree with many theoretical vertical-energy calculations of PDA. Singlet 
exciton relaxation studies on isolated polydiacetylene chains by subpicosecond pump-probe 
experiments have suggested that nonradiative singlet B^ exciton relaxation involves two Ag 



states in series it is well known that electron-phonon coupling is strong in PDA, 

we propose that lattice relaxation effects produce these Ag states. It is hoped that this will 
be resolved by studies of lattice relaxation explicitly within the Pariser-Parr-Pople-Peierls 
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Hamiltonian, which are now in progress. Similar studies by Barford on polyenes have found 
that the relaxation energies of the IBu and 2Ag can be reduced by as much as 0.3 and 1.0 
eV, respectively Thus, including lattice relaxation in the Pariser-Parr-Pople-Peierls 



model is expected to bring the 2Ag state below the IB^ in agreement with experiment. 
C Particle-hole Separation 

In order to help characterise the low- lying Ag and Bu states the particle-hole separation is 
calculated and gives an indication of the spatial extent of a given state: an explanation of 
how it is calculated is given in the Appendix. Suffice it here to say that we use the one- 
particle singlet excitation correlation function, which directly relates a hole in the valence 
band to an electron in the conduction band. In essence, if a particular state's particle-hole 
separation is seen to increase with system size, the electron and hole are unbound. However, 
if this quantity reaches a maximum, the Coulomb attraction between the electron and hole 
is strong enough to bind them together. In Fig. 8. the particle-hole separations are given 
as a function of system size in units of the average C-C bond length for the IBu, 2Ag and 
mAg states. It is clear that the IBu, 2Ag and mAg states are all bound excitons because 
the electron-hole separation reaches a maximum value with increasing chain length. These 



signify the composite particles mentioned in Section (|^). However, the particle- hole 
separation of the n5„ state increases linearly with system size: this suggests that the nB^ is 
either very weakly bound or unbound. We conclude therefore that for PDA molecules of up 
to 26 carbon atoms the electron-hole continuum (i.e., the band edge) is expected to be the 
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nBu, further strengthening the arguments of Section |^ (P). 

V CONCLUSIONS 

The DMRG calculation of a suitably parameterised Pariser-Parr-Pople Hamiltonian for a 
rigid geometry can be used to describe the electronic structure of short or long chain polydi- 
acetylene polymers. The role of electron-electron interactions is fundamental in the low-lying 
excitations, and this is borne out in our calculation. We have found and characterised the 
essential states: these are the lAg, IB^, mAg and n5„. Our IB^ excitation energies are 
within a few tenths of an eV of the experimental ones: for short chains, the calculated E 
(IBu) agrees with the experimental value, within solvation effects (~ 0.3 eV). The nBu is 
found to be above the charge gap, and hence it is continuum state; the 2Ag and mAg 
are not, and hence are bound excitons. For large chain lengths the nB^ tends towards the 
charge gap as expected, strongly suggesting that the nBu is the conduction band edge. We 
found our calculated conduction band edge to agree reasonably with the experimental value 
of ~ 3.0 eV quoted in the literature. Accounting for the strong polarisation effects of the 
medium and polaron formation gives our calculated E^{nBu) ~ 3.6 eV, with an exciton 
binding energy of ~ 1.0 eV. Our 2Ag is calculated to be above the li?„, which does not agree 
with relaxed-transition experiments: it is hoped that electron-lattice relaxation will adjust 
the position of the 2Ag to correct for this. Particle-hole separation data were used to help 
characterise the low-lying excitation in PDA. These data further suggest that the li?^, 2Ag 
and mAg are bound excitons, and that, conversely, the nBu is an unbound exciton. 
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This work is further evidence that the essential states mechanism is adequate for de- 
scribing the electronic properties of model vr-conjugated polymers. Future work will include, 
as already discussed, the incorporation of electron-lattice relaxation, and a more thorough 
treatment of the triple bond using a ZINDO Hamiltonian. 



We thank David Yaron for discussions. One of the authors (A. R.) is supported by the 
EPSRC (UK) . The calculations were performed on the DEC8400 at the Rutherford Appleton 
Laboratory. 



To measure the particle-hole separation we use the particle-hole correlation function, intro- 
duced in Refs. [48] and [49]. Underlying this approach is the assumption that a particle- hole 
pair corresponds to the promotion of an electron from the valence (or HOMO) band to the 
conduction (or LUMO) band. 

Without loss of generahty, the atomic basis can be transformed to a local molecular 
orbital basis. The local molecular orbitals correspond to the bonding (|1)) and anti-bonding 
1 2)) combinations of the atomic orbitals on each dimer, where a dimer is a double or triple 
bond. The local molecular orbitals delocalise via the single bonds. Then, the creation 
operator. 
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promotes an electron from |l))on dimer j. to |2)) on dimer i. However, since the local 
molecular orbitals are not exact Bloch transforms of the band molecular orbitalQ it is also 
necessary to include the hermitian conjugate of Sij. 

We now define the exciton correlation function with respect to the ground state, as: 

C!,{\n)) = {n\iSl + S,,)\l'Ap (A2) 

and the mean square of the particle-hole separation is: 

((^-jr) = ^^^|=j|^ (A3) 

In practice, we do not consider all combinations of i and j, but restrict ourselves to i and j 

symmetrically spaced around the middle dimer. 

^ The exact Bloch transforms of the band molecular orbitals arc Wannicr molecular orbitals. The Wannier 

orbitals of the valence band are predominately |1)), with a small admixture of |2)) from neighbouring dimers. 
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FIGURES 

FIG. 1. Chemical structure of a diacetylene trimer molecule, (a) Schematic representation of 
the orbitals giving rise to the different structures. The shaded orbitals give rise to cr-bonds, 
while the unshaded orbitals contain 7r-electrons in 2p orbitals. The orbitals all overlap 
along the length. The molecular plane is in the x-y plane, (b) Acetlyene structure, (c) 
Butatriene structure. 

FIG. 2. The idealised chemical structure of the series of PDAs used in this work. Mca 
represents three methyl structures. C„[A^] unambiguously represents the structures, where n 
is the number of carbon atoms and N is the number of monomer units. 

FIG. 3. Resonance integrals used for PDA with bond alteration in the Zth unit cell, ^^is 
the resonance integral for the single, tt for the triple, and td is for the double bond. The 
periodicity of polydiacctylcne is given by 4a, where a is the undimerized bond length. 

FIG. 4. Schematic energy band structure of polydiacetylene in the reduced Brillouin 
zone for the acetylene structure. E'^ is the energy gap within the valence bands. is the 
energy gap between the upper valence band and the conduction band. Numerical estimates 
are from the parameters of Table 2.0 (see text). 

FIG. 5. The four essential states. 

FIG. 6. Excitation energies for the 'iAg, mAg and n5„states for the PPP model as 
function of the inverse of the number of carbon atoms, 1/n. The li?„experimental excitation 
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energies are included for comparison. 

FIG. 7. The states contributing to the non-hnear properties of PDA and the important 
one-photon transitions between them. The dipole moments are shown for a chain of 26 
atoms. 

FIG. 8. Particle-hole separation for the IB^^ 2ylg, mAg and nB^ states in units of the 
average C-C distance as a function of the number of the number of carbon atoms, n. 
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d(C=C) 



D(C=C) 



D(C-C) 



sp angle' 



PDAs 

PDAs used in 
this paper 



PTS-6 
P-MBS 
P-THD 
Ci4[3] 

C22[5] 



d 
H 
H 



1.191 
1.195 
1.205 
1.202 
1.208 



1.356 
1.364 
1.426 
1.329 
1.360 



1.428 
1.424 
1.359 
1.427 
1.412 



121.9 

119.1 
123.0 
122.9 



a bond angles in degrees between an sp carbon atom and the next sp atom. 

b R=-CH2-OTos at 295 K. 

c R= -CHj -SOa-Ph-pOMe at 295 K. 

d R= -CH2-N(Ph)2 at 295 K. 
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a-bond 
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,c-c=c-c' 



H 



C.-C=C-C' 



R 



(b) 



R 



R C=C=C=C. 

^C=C=C=C'' 



C^C=C=C 



(c) 



Fig. 1. 



to 


2.539 eV 


a 


4.593 eV/A 


ro 


o 

1.4 A 


U 


10.06 eV 


ts 


2.4494 eV 


td 


2.7939 eV 


h 


3.4346 eV 
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Fig. 2. 



Chain 
Lenglh 


lAg - IBu 


IBli - mAg 


mAg - nBu 


n 




Dipole Moment 


AE 


Dipole Moment 


AE 


Dipole Moment 


AE 
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1 


1.5125 


5.1633 


2.2657 


8.8850 


2.5579 


14.0138 
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4.3072 
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7.5846 


2.7554 
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14 
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2.6615 


3.8860 


3.9430 


6.6115 


4.7893 


9.0628 


18 


4 


3.0926 


3.6483 


4.5245 


6.1595 


5.4395 


8.2927 
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5 


3.4679 


3.5015 


3.4514 


5.8547 


4.1472 


7.8628 


26 


6 


3.8018 


3.4049 


3.5629 


5.6994 


3.8258 


7.5910 
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Fig. 3. 
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(eV) 


IBu Theory 
Exp. 
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2.5113 


Charge gap, A 


5.7351 
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Figure 4 
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Fig. 7 
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